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Abstract. Single crystals of CeFeAsO, large enough to study the anisotropy of 
the magnetic properties, were grown by an optimized Sn-flux technique. The high 
quality of our single crystals is apparent from the highest residual resistivity ratio, 
RRR w 12, reported among undoped i?FeAsO compounds {R ~ rare earth) as well 
as sharp anomalies in resistivity, specific heat, C(T), and thermal expansion at the 
different phase transitions. The magnetic susceptibility x(^) presents a large easy- 
plane anisotropy consistent with the lowest crystal electric field doublet having a 
dominant Fg character. Curie- Weiss like susceptibilities for magnetic field parallel and 
perpendicular to the crystallographic c-axis do not reveal an influence of a staggered 
field on the Ce site induced by magnetic ordering of the Fe. Furthermore, the standard 
signatures for antiferromagnetic order of Ce at T^^ = 3.7 K observed in x(T) and C{T) 
are incompatible with a Zeeman splitting A « 10 K of the CEF ground state doublet 
at low temperature due to the Fe-magnetic order as previously proposed. Our results 
can be reconciled with the earlier observation by assuming a comparatively stronger 
effect of the Ce-Ce exchange leading to a reduction of this Zeeman splitting below 15 
K. 
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Since the discovery of superconductivity at = 26 K in fluorine doped LaFeAsO 
[1] and the increase of the superconducting transition temperature Tc up to 55 K in 
SmFeAsO [2] the layered FeAs systems and related compounds have been attracting 
considerable attention. In these compounds superconductivity appears close to the 
disappearance of an antiferromagnetic (AFM) state upon doping or application of 
pressure, suggesting an unconventional superconducting state in analogy to those 
observed in high-T^ or heavy-fermion superconductors. The AFM state below T^^ 
in the undoped compounds seems to be of itinerant nature and is connected with a 
large structural distortion [3l H]. While these unusual features are related to the 3d 
Fe-electrons, the CeFePnO [Pn = P, As) systems provide in addition strong electronic 
correlation effects due to the interaction between the Ce 4/-electrons and the conduction 
electrons. Thus, CeFePO was shown to be a paramagnetic heavy fermion metal, likely 
close to a ferromagnetic instability of the Ce 4/-system [5], while results on CeFeAsO 
indicate AFM ordering of Ce in a stable trivalent state [6l [7] . Recent results prove that 
the interaction between Fe-magnetism and 4/-magnetism in CeFeAsO differs strongly 
from that in the other i?FeAsO compounds {R = rare earth). Thus, Chi et al. [8] 
observed in inelastic neutron measurements on CeFeAsO a splitting of the Kramers 
doublets at the onset of the Fe-magnetism due to an internal field induced by the AFM 
Fe-order. This is surprising, since the high symmetry of the Ce site, respective to the 
columnar magnetic Fe structure, should result in the cancellation of Ce-Fe-Heisenberg- 
type exchange fields. /iSR results [9] evidence a larger hyperfine field at the muon 
site in CeFeAsO compared to all other i?FeAsO, which was attributed to an additional 
contribution originating from the staggered moment on the Ce site induced by the 
magnetic order of Fe. Very recently, investigations of the alloy system CeFeAsi-^P^O 
revealed a switching from an AFM to a ferromagnetic Ce-state at a; > 0.4 [10]. A likely 
origin for this unique behavior in CeFeAsO is a stronger hybridization between the /- 
states of Ce and the 3(i-states of Fe near the Fermi level connected with the inherent 
valence instability of Ce. 

Most investigations on the -RFeAsO systems have been performed on polycrystalline 
samples. Comparing the reported results with those obtained on 74Fe2As2 {A = Ca, 
Sr, Ba, Eu) samples, it is quite obvious that the quality of the i?FeAsO poly crystals 
is rather poor. Thus, the residual resistivity ratio (RRR) is usually less than 2 [7], 
compared with values extending up to, e.g., 32 in SrFe2As2 systems [11]. On the other 
hand, measurements on i?FeAsO single crystals are presently strongly limited by the 
very tiny size of the available single crystals, the largest length being typically less than 
0.3 mm [12]. This small size is a consequence of the most widely used preparation 
technique, based on a single crystal growth in NaCl flux (which has a low solubility for 
the involved elements) under high pressure, which puts strong limitations on the size of 
the crucible. 

We recently extended our Sn-fiux growth technique, which we developed for 
CeRuPO [I3], to the growth of CeFeAsO. After some optimization we obtained larger 
single crystals, which allowed the study of further properties, like e.g., the anisotropy of 
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the susceptibility. Further on, both a much higher RRR and narrower anomahes in the 
specific heat and in the resistivity at the magnetic and structural transitions indicated a 
much higher quality than the polycrystals reported so far. In this paper, we present the 
investigation and the analysis of the magnetic, thermodynamic, and transport properties 
of these CeFeAsO single crystals, focused on the magnetism of the Ce 4/-electrons. 

A typical example of a CeFeAsO single crystal is shown in Fig. [St, the edges are 
not very well defined due to the etching procedure after the growth. We found out 
that CeFeAsO is attacked by diluted HCl on a much higher rate than the CeTPO- 
compounds (T = Fe, Ru, Os, Co); therefore, we centrifuged the compound in the 
Sn-matrix before etching for only 10 minutes. The single crystals were characterized 
with energy dispersive X-ray analysis revealing a stoichiometric Ce:Fe:As content and 
confirming the presence of oxygen in the compound; furthermore, no foreign phases 
could be detected. We specifically looked at Sn incorporation, but did not observe any 
Sn line in the microprobe spectra of several 1111 single crystals. In the Sn-fiux crystal 
growth of AFe2As2 compounds, Sn-incorporation was observed to be quite severe for A 
= Ba, but much weaker for A = Sr, Ca [H]. 

Several X-ray powder-diffraction patterns on crushed single crystals recorded on 
a Stoe diffractometer in transmission mode using monochromated Cu-Kq, radiation 
(A = 1.5406 A) confirmed the PA/nmm structure type and the formation of single 
phase CeFeAsO. The lattice parameters a = 4.002(1) A and c = 8.647(2) A, refined 
by simple least-squares fitting, were found to be in good agreement with the reported 
structure data [U El [3, [15] . Susceptibility xiT) and magnetization M{H) measurements 
were performed in a commercial Quantum Design (QD) magnetic property measurement 
system (MPMS). Resistivity p(T) was determined down to 0.4 K using a standard four- 
probe geometry in a QD physical property measurement system (PPMS). The PPMS 
was also used to measure the specific heat C{T) with a standard heat-pulse relaxation 
technique. High-resolution thermal expansion measurements were carried out using 
a capacitive dilatometer built after [16], enabling to resolve relative length changes 
A/// > 10-10. 

The quality of our single crystals is demonstrated by the resistivity data presented 
in Fig. [IK. In this measurement the current was applied along the basal plane. The 
room-temperature value, psooK = 1.54mr2cm, indicates bad metallicity, followed by an 
upturn towards the structural transition at Tj*^ = 149 K. Below T|^° and T]^°, p{T) 
decreases drastically in a metallic way reaching p2K = 140 fiQcm at T = 2 K, giving 
the highest RRR ~ 12 among the -RFeAsO samples studied so far. T|^^ and T^'^ were 
determined from the derivative of the electrical resistivity. We are presently performing 
an analysis of the sample dependence of these transitions, which shall be presented 
in a forthcoming publication. For our best single crystal we found T|^® = 149 K and 
T^^ = 145 K. 

The high quality of our single crystals is further highlighted by the sharpness of 
the transitions at Tg^^ and T^^. The large size of the anomaly in p{T), a{T), and 
C(T) (not shown) confirms the presence of the AFM ordering of Fe in our single 
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Figure 1. (a) Resistivity as function of temperature for a CeFeAsO single crystal. 
The sharp drop at T|''',T^° = 150 K hallmarks the structural phase transition and the 
AFM ordering of Fe. (b) p{T) around the anomaly at T^'^. (c) The magnetoresistance 
at T = 2 K is linear up to 10 T and quite large. 

crystals. Furthermore, the pronounced drop in p(T) we observe at T^'^ and 
supports the absence of any significant Sn-incorp oration in our single crystals, because 
Sn incorporation in Sn-fiux-grown BaFe2As2 single crystals results in an increase in 
resistivity when cooling below T^*^, Tg'" in contrast to a decrease for samples grown 
from self-flux (see e.g. [H]) and for polycrystalline samples; also recent results show an 
increase in resistivity when cooling below Tj^*^, T^'^ in some large LaFeAsO single crystals 
grown from NaAs-flux, this increase being attributed to defects [I7j. In Fig. [T]d, p(T) 
is enlarged for T < 10 K, where a distinct anomaly is visible at the AFM ordering 
temperature of Ce, T^^ . At this anomaly, p{T) tends to slightly increase before 
decreasing with lower T. The former suggest the opening of a small gap at the Fermi 
energy, the latter is because of reduced spin-spin scattering in the AFM ordered state 
and typical for Ce-systems. Preliminary measurements of the magneto-resistance p{H) 
for magnetic field applied along the c direction, shown in Fig. [T]^, reveal a huge positive 
magneto-resistance, reaching Pfj,oH=iii / PtJ,oH=o = 4.3. For fields up to poH = 10 T, 
p{H) is linear followed by a steeper increase at higher fields. Further measurements 
with different magnetic field directions are necessary to clarify the origin of this greatly 
enhanced magnetoresistance and its peculiar field dependence in CeFeAsO. 

The magnetic susceptibility of CeFeAsO is dominated by the contribution of the 
well localized, trivalent Ce. While the paramagnetic susceptibility of Fe-systems is 
generally very isotropic due to a small spin orbit coupling (SOC), the large SOC of 
rare earth elements can lead to a huge anisotropy of the magnetic properties. That is 




Figure 2. Anisotropic susceptibility of CeFeAsO for _ff ± c and H \\ c. The AFM 
order of the Ce-ions at T^^ = 3.7 K is clearly visible for both crystallographic directions. 
Inset: Inverse susceptibility as function of temperature measured in a magnetic field 
of iiqH = 5T. The open symbols (black) represent the measured data and the closed 
symbols (red) denotes the 4/-part of the inverse susceptibility obtained by subtracting 
XLaFeAsO reported in Ref. [H]. A Curie- Weiss fit above 150 K gives fi^^ = 2.45 for 
H ± c and i^lg = 2.6 for H \\ c (blue fines). 



the case for CeFeAsO as evidenced by our susceptibility data on a single crystal. In 
Fig. [2] we compare x{T) for H A. c and H \\ c. x± increases continuously below 300 
K and reaches a value of 1.16-10~^m^/mol at the maximum at T^^ . In contrast, x\\ 
is almost constant between 300 and 100 K. Although it starts to increase significantly 
below 100 K, it reaches only a value of 0.2-10~^m^/mol at T^'^ , a factor of 6 smaller than 
within the basal plane. The Curie- Weiss plot x{T)~^ versus T allows for a more detailed 
analysis (inset of Fig. [2]). While xl^ shows an almost straight line from T^"^up to the 
highest temperatures with an effective moment ii^q^'^^^^ = 2.71 ^b, slightly larger than 
that expected for a free Ce^"*" ion, increases linearly only up to 80 K, and then stays 
almost constant. For a more precise determination of the 4/-contribution we subtracted 
the reported susceptibility of LaFeAsO [18], XLaFeAsO, since the Fe and valence electron 
contributions are not expected to differ significantly between the La- and the Ce-based 
compounds. Fits (blue lines) to these xl/ = (xceFeAsO — XLaFeAso)"^ data (closed red 
symbols in the inset of Fig. [2]) give an effective moment of /i^ = 2.45 /x^ and a Weiss 
temperature of 0^ = IIK for if ± c (fit above 150 K), while for H \\ c we get 
/^eff-highT — '^■Q f^B and Oiligijj, = — 580K for the fit above 150 K but /i|!g_io„r — 0.96 /is 
and Qiovjt — — 3.6K for the fit below 70 K. 

The anisotropy of the susceptibility in rare earth compounds is governed by the 
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crystal electric field (CEF) and to a less extent by the anisotropy of the exchange 
interactions respective to the different components of the local moment. The strong 
anisotropy, we observed at high temperatures, point to a huge CEF anisotropy in 
CeFeAsO; thus, we can in a first approach neglect the exchange anisotropy. Then, the 
leading CEF parameter 5° (Stevens coefficient) can be estimated from the difference 
between 0^ and By, determined at high temperatures ([T9| 120]): 

^°-<^--^') 3.(2J-l)°.(2J + 3) « 

In our case the corrections for the Fe-contribution in x\\ l^ad to a very large 
uncertainty in OhjghT) ^ind therefore the value 5° = 62 K obtained from our data 
only corresponds to an order of magnitude estimate. In this sense it is in reasonable 
agreement with the value i?2 = 37K proposed by Chi et al. [8] on the basis of neutron 
scattering results. The very well defined slope of versus T below 70 K allows an 
estimate of the c-axis effective moment of the ground state doublet, and thus, a further 
insight into the CEF scheme. For Ce"^"*" in a local tetragonal surrounding with a large, 
positive, and dominant B2 CEF parameter, the CEF ground state is expected to be the 
Fg doublet (| ±1/2)). This was indeed proposed to be the case in CeRuPO [13] as well 
as in CeFeAsO above the ordering temperature of Fe [8]. Fg bares a strongly reduced 
c-axis saturation moment of fi^^^ = 0.43 fiB and a three times larger a;|/-saturation 
moment fi^^^ = 1.29 ^b- Accordingly, one expects at lower temperatures a strongly 
reduced effective moment along the c direction, of only = 0.74 /i^, while in the 
basal plane the effective moment = 2.23 yUe is only slightly reduced compared to 
the value for the full Ce^"*" state. Our experimental results below 70 K confirm such 
a strong reduction along the c direction and almost no reduction for the basal plane, 
and therefore strongly support that the CEF ground state doublet has a dominant Fg 
character. The orthorhombic distortion at = 150 K shall induce a slight mixing of 
the I ±|) and | ±|) wave functions in the Fg state, which might be the origin of the 
slight difference between the Fg and the observed fif^. In the tetragonal phase above 
150 K, the excited CEF levels corresponding to a Fg ground state are two doublets, 
F^^^ = cos(a) I ±3/2) ± sin(a) | ^5/2) and Pf ^ = -sin(a) | ±3/2) ± cos(a) | ±5/2). Chi 
etal. proposed a rather large (negative) mixing angle a = —25.6° [8]. We calculated the 
susceptibility for such a CEF scheme assuming the excited levels to be at the energies 
observed in neutron scattering experiments but allowing for different mixing angles a 
(only the absolute value of a is relevant for x(T)). We found that x\\ is quite sensitive 
to |a|: While the absolute value of 1/X|| at 300 K is almost independent of |a| and 
matches the experimental value nicely, the slope, d(l/x||)/dT, at 300 K decreases with 
increasing |a| and becomes negative for |q;|>20°. Since we observe a positive slope, 
our susceptibility data suggest a smaller value of |a| than proposed on the basis of the 
neutron results, but unfortunately the large uncertainty in the experimental x\\ at high 
temperatures does not allow for a precise determination. 

Inelastic neutron scattering [8] and /iSR data [9] indicate a splitting of the CEF ground 
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Figure 3. (a) Anisotropic magnetization of CeFeAsO at T = 2K. (b) T"^^ determined 
from x{T) measurements (Fig. [2|) as function of the external magnetic field, (c) 
Photograph of a CeFeAsO single crystal, the scale labels 1 mm. 



state doublet and a corresponding Ce polarization, attributed to an Fe-induced internal 
field at the Ce site below T^°. Our susceptibility data do not reveal any effect of this 
splitting and polarization, neither in the T- nor in the if- dependence, which is at first 
sight surprising. We shall discuss this point in details later on. 

Below T^^ = 3.7 K, both xj_(T)and X\\{T) show a pronounced decrease, indicating 
AFM ordering of Ce. Here again, no evidence for an Fe-induced polarization of Ce 
is obvious; the x{T) dependence looks like in any standard Ce^^ system ordering 
antiferromagnetically. With increasing applied field, Tj^-^shifts to lower temperatures, 
more strongly for field in the easy plane than for field along the hard c direction, as 
expected for a classical antiferromagnet. 

A preliminary sketch of the H — T phase diagram is shown in the inset of Fig.[3}3. 
In the main part of this figure, Fig.[3h., we present the magnetization of a CeFeAsO 
single crystal at T = 2K, i.e., in the AFM ordered state, for both crystallographic 
directions. For H c, M{H) increases linearly for fields fioH < 2 T, followed by a 
small metamagnetic transition around fioH = 2T. The origin of this metamagnetic 
transition is presently not settled and thus deserves further magnetization as well as 
microscopic measurements. For fioH > 2T the magnetization increases again linearly, 
reaching an absolute value M = O.S/i^/Ce at fioH = 5T, which is about 60% of the 
basal plane saturation moment, /x^^* expected for a Tq doublet. Assuming a linear 
extrapolation of M{H) this saturation value would be reached at /io-f^sat ~ 7.8 T. For 
H II c, M[H) presents a perfectly linear behavior up to fioH = 5T with a much 
smaller slope dM/dB = 0.019 fJ,B/T- These results are in agreement with the easy- 
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plane anisotropy of CeFeAsO determined by the susceptibility measurements discussed 
above. 

Within a simple molecular field model, defining the molecular field B^y induced by the 

M as 

xy 



antiferromagnetic Ce-Ce coupling constant A;^™ as 



^xy ~ ^xy ' '^xy i (2) 

with m^y being the ordered Ce moment (in analogy to equation 12 in |[9]), one can 
estimate A^™ either from the initial slope dM/dB = 0.11 i^b/T 

A^™ = l/{2dM/dB) = 4.5 T//iB (3) 

or from the saturation field i?sat 

= ^sat/(2/.r;) = 3.0 T//iB. (4) 
Within the same approximation, an independent estimation can be obtained from T^^ , 

= ikn/fiByT^^/if^f^r = 3.4 T//iB, (5) 

with fi^^y in units of /ib. All these values are in reasonable agreement, suggesting a 
rather conventional type of magnetic ordering. Accordingly, the initial slope dM/dB 
for field along the c-axis leads to \^^^ = 26 T//iB- The fact that Ce orders in the basal 
plane and not along the c direction gives an upper limit for A^^^ 

Af ^ < (kn/fiB) ■ T'J/if^rr = 30 T//iB. (6) 

Comparison with X^™ suggests a huge exchange anisotropy, with a much larger 
exchange for the z component. However, this merely refiects the anisotropy of the g 
factor, since in the molecular field approximation A is proportional to g~^, (e.g. |21j). 
which is a factor 9 higher for the z direction compared to the x|/-direction. Thus, the 
product Xi- gf, which is proportional to the exchange parameter in a Heisenberg picture, 
is slightly lower for the z component than for the basal plane component. 
Next, we focus on the specific heat measurements around T^^ , plotted as C^^ /T vs T 
in Fig. H^. The 4/-contribution, C"^-^, to the specific heat was obtained by subtracting 
the reported specific heat data of LaFeAsO [22] which, however, is small compared to 
the large anomaly due to the onset of AFM ordering of the Ce-ions. This sharp A- 
type anomaly refiects the high quality of our single crystal compared to measurements 
on polycrystals [6]. The value C^ax = 12.6J/(mol K) at the peak at Tj^'^is almost 
identical to the value 12.5 J/(mol K) calculated in a simple mean field model for a 
doublet system. From the anomaly in C^^/T the AFM transition temperature can be 
accurately determined yielding = 3.7 K. The anomaly shifts to lower temperatures in 
an applied magnetic field, similar to what is observed in the susceptibility measurements 
(see Fig. [2] and [3b). In zero field, C^^ /T decreases quadratically with temperature below 
T^^ , supporting that the excited magnons are of AFM nature. Below T = 0.8 K, C^^ /T 
levels off at a constant Sommerfeld-coefficient, 70 = 50 mJ(mol K^)~^, which is enhanced 
by more than one order of magnitude compared to LaFeAsO (79 = 3.7mJ(mol K^)~^ 
[22]). This enhancement can therefore be attributed to correlations of the 4/-electrons. 
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Figure 4. (a) 4/-part of the specific lieat of CeFeAsO around . Tlic anomaly 
shifts to lower T in an external applied magnetic field, (b) Linear thermal expansion 
coefficient ai as a^/T for single crystalline CeFeAsO measured along an unspecified 
in-plane {ah) direction, a_\_, and along the c-axis, a\\. 



Whether or not these correlations are also entangled to the 3(i-electrons can presently 
not be answered and deserves further more microscopic studies. 

The magnetic entropy was obtained by integrating C^^ /T over T and is presented 
in Fig. [5b- The steep increase at T^^ together with absolute values of 80% of i?ln2 at 
T = 8 K confirm that most of the 4/-degrees of freedom of the lowest CEF doublet are 
involved in the AFM order. 

Finally, in Fig.Hb we show the results of the uniaxial thermal expansion coefficient 
a(T) = l~^{dl/dT) along a non-specified direction within the a6-plane, a± and along 
the c-axis, a\\, for CeFeAsO below 8K in a representation a/T vs. T. The data are 
dominated by pronounced and strongly anisotropic phase transition anomalies at the 
AFM ordering of the 4/-moments at T^/ = (3.7 ± 0.1) K. Both the sharpness of the 
transition and the pronounced anisotropy confirm the high quality of the single crystal 
studied here. As shown in Fig.Hb, the negative discontinuity in a\\ is preceded by 
an anomalous negative expansivity, unexpected for a lattice contribution, but likely 
to originate from the Ce-4/ degrees of freedom. Its pecuhar T-dependence suggests 
a relation to the decrease of the Zeemann splitting discussed later on. The c-axis 
behavior contrasts with a small and positive a± revealed above T^^ . Note that a± can 
be significantly affected by the domain structure which forms upon cooling through the 
tetragonal to orthorhombic phase transition at Tj" = 150 K. Details of the expansivity 
results at T^^ will be published elsewhere. 

Combining thermal expansion with specific heat data (Fig.Hj) and making use of 
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T{K) T{K) 

Figure 5. (a) Calculated (see text) and measured magnetic susceptibility for 
temperatures below 30 K. (b) Magnetic entropy calculated by integrating C^^ /T over 
T. (c) Specific heat plotted as C'^^ vs T; for comparison the specific heat calculated 
from the proposed CEF scheme [5] is drawn as a line (see text). 



the Ehrenfest relation. 

dTtI 



Vmoi ■ TtJ ■ Aa,/AC, (7) 

an estimate of the pressure dependence of T^^ can be made. Using the molar 
volume Vmoi = (4.17 ■ 10"^) mVmol, Aa\\ = -(10.5 ± 0.5) ■ 10"^ K"^ estimated from 
Fig.Sb and AC^/ = 12.6 Jmol-^K-\ one finds dT^^/dP\\ = -0.013 K/kbar. Although 
a quantitative estimate of the in-plane pressure dependence is precluded by potential 
domain structure effects, the large positive discontinuity in a±_ at T^^ in Fig.HjD indicates 
a positive in-plane pressure effect on T^^ . Since a very small pressure dependence of 
T^^ is expected under hydrostatic pressure conditions, (substituting large As by smaller 
P left T^'^ almost unaffected up to x < 0.3 in CeFeAsi_j;Pa;0 [10]), a significant in-plane 
anisotropy in dT'^^ / dPa^b can be inferred. 

According to the results and discussion presented so far, CeFeAsO behaves like a 
standard trivalent Ce^"'' system ordering antiferromagnetically at T^^ = 3.7 K, with a 
strong easy-plane CEF anisotropy. There is neither evidence for a large Kondo effect 
as in CeFePO, nor direct evidence for an influence of the SDW ordering of Fe. More 
specifically, the observation of a well-defined AFM transition at Tj^'^ associated with a 
large entropy seems to be in strong contradiction with a Zeeman splitting of the CEF 
ground state of 0.93 meV = 10.8 K as proposed by Chiet al. [8], on the basis of inelastic 
neutron scattering measurements. On the one hand, for such a splitting the polarization 
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of the ground state already amounts to 92 % at 4.5 K, i.e., the system would already be 
almost completely polarized at a temperature 20% above the observed T^-^ . Thus, there 
would not be much degrees of freedom left for a magnetic transition. Accordingly, most 
of the entropy of i?ln2 associated with the CEF ground state doublet would have been 
released at higher temperatures, the entropy at 4 K left for a transition would amount 
to only 34 % of Mn2, far below the observed value 55 % of Mn2. Furthermore, the 
4/ specific heat should be quite different from that observed experimentally. We have 
calculated the specific heat using the appropriate Schottky-model for 4 non-degenerated 
energy levels at = 0, 0.93, 16.9, 20.1 meV as proposed by Chi et al. [8]. In Fig. [5t 
the resulting specific heat curve is compared to our experimental data as C^^ vs. T 
around T^^ . The calculated curve shows the expected broad Schottky-type anomaly 
with maximal values Cmax = 3.65J/molK at T = 4.5 K. This is in strong contrast to 
the A-type anomaly with absolute values C^-^ = 12.6J/molK at T^^ and the almost 
T-independent C(T) above 5 K. 

In order to shed more light into this apparent contradiction, we performed a more 
profound analysis of our results. The first point we looked at in more detail was how 
the Ce susceptibility is affected by an internal field due to the SDW ordering of Fe. We 
assumed a Fg CEF ground state doublet, which corresponds to an effective S = j state 
with an anisotropic g factor g^y = 2.57 and gz = 0.86, under the influence of both an 
internal field i?int and an external field i?cxt- Maeter et al [H], give strong evidence that 
the internal field on the Ce site due to the Fe ordering is a staggered field parallel to the 
c-axis. We followed this proposition and took 5™* = 18.75 T as required for a splitting 

A = B'f -gz-fiB/kB = 10.8 K. (8) 

We assumed -B"^* to be T-independent since we were interested in the low-temperature 
behavior (T < 30 K), where the Fe-moment is already saturated. For an external 
field along the c direction, the magnetization and the susceptibility can then easily 
be calculated using the usual Brillouin function Bi/2{x) for two different Ce-sites, one 
with Sext and Sjnt parallel and one with _Bext and Bi-^t antiparallel. We found that B^^t 
has almost no effect above 15 K, but below 15 K it starts to reduce x\\y such that the 
c-axis susceptibility should show a maximum at T = 7 K followed by a pronounced 
drop below 6K (Fig.[5K, upper part). This can easily be understood from the Brillouin 
function since the susceptibility is now related to the slope at 

X = ±Ant ■ gz ■ f^B/{2kBT) (9) 

instead of the slope at the origin x = 0. As long as x is in the linear part of the 
Brillouin function, there is no difference, but once one reaches the saturation region the 
slope vanishes. We also studied x±] for a Zeeman splitted doublet due to Bi^t along the 
c-axis, x± corresponds to a Van Vleck susceptibility due to mixing between the | +1/2) 
and the | —1/2) states: 



Rare earth magnetism in CeFeAsO: A single crystal study 



12 



This leads to a constant x± low T, but as soon as the temperature gets closer 
to A/ks, the excited state gets populated and x± merges in a Curie law with the 
usual moment (Fig.[5|L, lower part). The calculation revealed that the departure from 
the Curie law is only significant below 5 K, thus in a temperature range where the 
experimental data are affected by the antiferromagnetic ordering. In conclusion, this 
analysis clarifies the absence of any effect of i?int on x±{T) above T^'^and on x\\ (T) above 
15 K, but confirms that the observed continuous increase of x\\ {T) below 10 K down to 
T^^ is in clear disagreement with the suggested splitting, which should lead to a visible 
maximum. This cannot be repaired by changing the orientation of Ant; since then an 
equivalent maximum would appear in x±(T). Therefore, both the specific heat and x\\ 
contradict a splitting of the CEF ground state of the order of 11 K below T < lOK. 

Thus, we are left with a contradiction between the clear evidence for a splitting at 
higher temperatures in the neutron and /xSR results and the absence of the expected 
effects of this splitting in x\\ and C{T) below 10 K. This contradiction can be resolved by 
assuming a reduction of the splitting towards lower temperatures. A closer inspection of 
the neutron data of Chi et al. [8] supports such an assumption. The low-energy spectra 
in the inset of Fig. 2a of this reference reveal a quite sharp excitation at 1.05 meV (= 
12.2 K) for T = 30K and T = 96 K, but for T = 7.8 K the peak has shifted to 0.85 
meV (= 9.9 K), and seems to have broadened on the low-energy side. In order to get a 
further idea about the possible magnitude of a T-dependent splitting, we assumed that 
above 5 K the system can be fully described by two singlets separated by a T-dependent 
splitting A(T), and calculated A(T) from the measured entropy S^^ , 

— A(T) 

S"^^ = R{pQ\npQ + pilxipi) with Pi/po = exp{ }. (11) 

The results are plotted in Fig.E] together with the peak positions in the low-energy 
inelastic neutron data. Both data sets consistently indicate a decrease of the splitting, 
especially below 10 K. A T dependent internal field was already introduced by Maeter 
et al. In order to get a better fit to their /iSR data, they considered a further 
exchange field on the Ce-site due to the polarization of the Ce itself: 

where the exchange field due to the Fe was assumed to be T-independent below 30 K, 
^Fc-Cc _ j^^^ j^^ corresponds to a Ce-Ce coupling constant for a specific symmetry of 
the Ce exchange, and m^^ is the polarization of the Ce-site. 

Maeter et al. [9] obtained the best agreement with the yuSR data for Bq = 26.8 T 
and Jsz = — 24.3T//iB. Interestingly, this coupling parameter is quite close to the cou- 
pling constant we estimated from x\\ below Tn (see discussion of susceptibility data). 
The negative sign of J^^ leads to a reduction of the total internal field with decreasing 
temperature, because of the increasing polarization of the Ce in the Fe exchange field. 
Accordingly, the splitting of the Fg ground state becomes T dependent. 
We calculated A(T) using the parameter of Maeter et al. [9] (m^^ was calculated in an it- 
erative process using the Brillouin function), and plotted the results in Fig.[6l where they 
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Figure 6. Gap energy A as a function of temperature estimated from neutron 
scattering [8] and entropy data. The open triangles are the values calculated with the 
model and parameters of Maeter et al. [9]. In a modified calculation (closed triangles) 
the Ce-Fe exchange field was slightly reduced to reconcile experimental and calculated 
values. 



can be compared with the values deduced from neutron and specific heat experiments 
as described above. While the calculated T-dependence shows the expected decrease 
of A(T) with decreasing T, the absolute values are significantly above the values we 
deduced from experiments. A nice agreement can be obtained by just reducing the Fe- 
Ce-exchange field to Bq = 23 T. The effect of the Ce-Ce exchange is then comparatively 
stronger, and is already quite sizeable at 10 K. The broadening on the low-energy side of 
the excitation in the low-energy neutron spectra at T = 7.8 K, as well as further analysis 
of the T dependence of the specific heat and of the susceptibility indicate that the effect 
of the Ce-Ce exchange increases more strongly towards lower T than suggested by this 
simple model. One would e.g. expect that the inherent tendency of Ce towards an own 
magnetic ordered state leads to short-range correlations above the ordering tempera- 
ture, which would be associated with a fluctuating transverse fleld, not accounted for 
in the simple 'static' model. In summary, our analysis of the susceptibility and speciflc 
heat data together with a reassessment of the neutron data indicate a strong effect of 
the Ce-Ce exchange at lower temperatures, competing with the exchange fleld induced 
by the Fe-ordering, and leading to a decrease of the total internal fleld at the Ce site 
starting already below around 15 K. 

In conclusion, we have grown larger single crystals of CeFeAsO using a Sn-flux tech- 
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nique and investigated the magnetic properties related to Ce by means of x(T), M{H), 
C{T,H), p{T,H), and <y{T) measurements. In contrast to earlier results on polycrys- 
tals, the basal plane resistivity in our single crystals not only shows a sharp drop at 
the structural (and the antiferromagnetic) transition at T = 150 K, but continues to 
decrease down to the lowest temperature, resulting in a high RRR ^ 12. This ob- 
servations evidence on the one hand the presence of a well-defined metallic state in 
CeFeAsO at low temperatures, and on the other hand the high quality of our single 
crystals, which is also demonstrated by the sharpness of the anomalies in the resistivity, 
specific heat, and thermal expansion at the different transitions. The magnetic suscepti- 
bility, which is dominated by the contribution of the well localized 4/ electrons of Ce^"'', 
presents a strong easy-plane anisotropy, with predominantly Fg (| ±1/2)) character. At 
T^-^ = 3.7 K, maxima in the susceptibility along both directions as well as a sharp A-type 
peak in the specific heat and in the thermal expansion evidence a transition to an anti- 
ferromagnetic state. The large amount of entropy released at the transition implies that 
most of the degrees of freedom of the CEF ground state are involved in the formation 
of the antiferromagnetic state. This is also confirmed by the peak value of C (T) being 
very close to that expected for a mean-field transition. This large amount of entropy 
at Tj^'^as well as the continuous increase of x\\ below 15 K down to Tj^-'^are incompatible 
with a T-independent Zeeman splitting of the CEF ground state A = 10.8 K as sug- 
gested from neutron measurements. The analysis of the specific heat and reassessment 
of the neutron data indicate a significant reduction of the splitting below 15 K. This 
reduction can certainly be attributed to the growing effect of Ce-Ce exchange, which 
is competing with the Fe-Ce exchange and thus leads to a reduction of the molecular- 
field at the Ce-site. Although this is likely a precursor effect of the Ce-ordering at 
T^^ = 3.7 K, it starts to become visible at a much higher temperature, around 15 K. 
The 4/-electrons in the ordered state show pronounced correlations, reflected in an en- 
hanced 7o = 50mJ/molK^, more than one order of magnitude higher than in the closely 
related compound LaFeAsO without 4/-electrons. The observation of a heavy fermion 
state connected with strong ferromagnetic Ce-Ce exchange in CeFePO suggests a very 
complex interplay between 4/ and 3d physics in the alloy CeFeAsi-ajP^jO, both on the 
level of exchange effects and on the level of strong correlation effects. This certainly 
deserves further investigations. 
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